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ABSTRACT 

We have detected the T Tauri star, DO Tauri, in a 0.6"-resolution VLA map of 
43.3 GHz (A = 7 mm) continuum emission. The 43 GHz flux density lies on the 
same power-law slope defined by 89 to 232 GHz measurements, F u oc v a with index 
a = 2.39±0.23, confirming that the 43.3 GHz emission is thermal radiation from 
circumstellar dust. Upper limits to the flux densities at 8.4 and 22.5 GHz constrain 
the contribution of free- free emission from a compact ionized wind to less than 49%. 
The dust emissivity index, /3, is 0.39±0.23, if the emission is optically thin. Fitting a 
model of a thin circumstellar disk to the observed spectral energy distribution gives 
(3 = 0.6 ± 0.3, consistent with the power-law derivation. Both values are substantially 
lower than is generally accepted for the interstellar medium, suggesting grain growth. 
Given the youth of DO Tau and the early evolutionary state of its circumstellar disk, 
this result implies that mm-size grains have already formed by the early T- Tauri phase. 



Subject headings: stars: individual: DO Tauri — stars: formation — circumstellar 
matter — planetary systems 
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1. Introduction 

At least 50% of T Tauri stars (TTs) appear to be surrounded by circumstellar dust disks 
(Strom et al. 1989; Beckwith et al. 1990, hereafter BSCG; Andre & Montmerle 1994; Henning & 
Thamm 1994; Osterloh & Beckwith 1995). Global disk properties can be inferred from models 
of spectral energy distributions (SED's) from infrared to millimeter wavelengths (Adams, Lada, 
& Shu 1987; Beckwith & Sargent 1993; Mannings & Emerson 1994). Masses and sizes are similar 
to those assumed for the early solar nebula, suggesting that the disks may be protoplanetary (cf. 
BSCG; Beckwith & Sargent 1993). However, the SED models rely on assumptions about disk 
morphology and radial structure, and about the nature of the constituent dust grains. 

Grain size and composition in these potentially planet-forming disks can be inferred from 
the spectral index, f3, of the dust opacity (cf. Pollack et al. 1994) if thermal radiation from 
grains in the disk is optically thin. Sub-arcsecond resolution is necessary to image disks directly 
and measure properties on spatial scales of <100 AU at the distance of the nearest star-forming 
regions. At wavelengths longer than 3mm, the emission is well into the Rayleigh- Jeans part of the 
Planck curve and very likely to be optically thin. Dust continuum radiation has been detected 
from a number of TTs at A = 3 mm (Sargent & Beckwith 1993 and references therein), but no 
thermal emission has been detected unambiguously at longer wavelengths (Mundy et al. 1993). 
The required spatial resolution and mJy sensitivities can now be achieved using the VLA at 
wavelengths of 7 mm. 

DO Tauri is a young TTs in the Taurus star formation complex at a distance of 140 pc (Elias 
1978; Kenyon, Dobrzycka, & Hartmann 1994). Estimates of its age and mass range from 1.6 to 
6.0 xlO 5 yrs and 0.3 to 0.7 M Q (BSCG; Hartigan, Edwards & Ghandour 1995), depending on 
the theoretical tracks used to place it on the H-R diagram. The spectral energy distribution is 
consistent with the presence of a ~0.01 Mq circumstellar disk (BSCG; Beckwith & Sargent 1991, 
hereafter BS; Mannings &; Emerson 1994). Asymmetric, blue-shifted, [OI] and [SII] forbidden 
line emission (Appenzeller, Jankovics &; Ostreicher 1984; Edwards et al. 1987; Edwards, Ray &: 
Mundt 1993) is resolved as an optical jet at PA 70° (Hirth et al. 1994). The jet is approximately 
orthogonal to the direction of linear optical polarization, PA ~ 170° (Bastien 1982), and to the 
long axis of CO (2— >1) emission detected in aperture synthesis images of DO Tau, PA ~ 160° 
(Koerner &: Sargent 1995). Kinematic models of the molecular line emission are consistent with 
the presence of a circumstellar disk that is centrifugally supported within a radius of 350 AU from 
DO Tau (Koerner 1994). 

Here, we report on sub-arcsecond images of the A = 7 mm emission from DO Tau which were 
made using the recently upgraded Very Large Array (VLA) of the National Radio Astronomy 
Observatory^. We have supplemented these measurements with continuum observations of DO 
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Tau at other wavelengths to sample the spectral distribution of emission from A = 1.3 mm to 
3.6 cm and improve our understanding of grain properties in the circumstellar material. 



2. Observations and Results 

The VLA was used to observe DO Tau in radio continuum emission at 43.3, 22.5, and 8.4 
GHz (A = 7 mm, 1.3, and 3.6 cm). The phase center was offset 1" in both RA and Dec from the 
stellar position of DO Tau (Herbig & Bell 1988), and the total bandwidth was 100 MHz in right 
and left circular polarizations. As for all observations discussed below, molecular line emission is 
negligible within the narrow band. Observations at 43.3 GHz were carried out on 1994 April 3-4 
with the inner seven antennas of the high-resolution A configuration, and on 1994 August 20 with 
10 inner antennas of the B configuration. Baselines up to 5.6 km provided UV coverage in the 
range 30-800 kA. On both dates, DO Tau was observed at 22.5 and 8.4 GHz using the remainder 
of the VLA's 27 antennas. UV coverage was 50-2700 kA at 22.5 GHz and 20-1000kA at 8.4 GHz. 
Absolute flux densities were calibrated using 3C48 and 3C286 with an estimated uncertainty of 
20%. At 43.3 GHz, gain calibration was accomplished with periodic observations of 0333+321 
with a measured flux density of 0.98 ± 0.06 Jy. At 22.5 GHz and 8.4 GHz, the gain calibrator was 
0400+258 with flux densities 0.65 ± 0.03 Jy and 0.83 ± 0.01 Jy, respectively. 

Data calibration and mapping used standard routines in the NRAO AIPS software package. 
Daytime atmospheric phase fluctuations during A array observations necessitated extensive 
editing and application of a Gaussian taper to the UV data, resulting in a 0.68" x 0.53" (FWHM) 
synthesized beam at PA —78°. This corresponds to 95x74 AU at DO Tau. Fig. la displays the 
CLEANed image of DO Tau at 43.3 GHz. An unresolved source with flux density 1.80 ±0.71 mJy 
is detected at the stellar position, a(1950) = 04 fe 35 m 24.19 s , 5(1950) = 26°04'54.5". The ±0.71 
mJy uncertainty includes rms variations in the map (± 0.35 mJy bm -1 ) and a possible 20% error 
in absolute flux calibration. At 22.5 and 8.6 GHz, DO Tau was not detected within the area 
encompassed by the 43.3 GHz synthesized beam to 3a levels of 0.76 and 0.17 mJy, respectively. 

Observations were made with the Owens Valley millimeter array at 89.2, 111.2, 221.5, and 
232.0 GHz (corresponding to A = 3.4, 2.7, 1.4, and 1.3 mm) between 1993 September and 1995 
March. Measurements at 89 GHz were made with six telescopes; four were used at 110 GHz, 
and five at 220 and 230 GHz. Overall UV-ranges were 5-60 kA (89 GHz), 5-25 kA (110 GHz), 
and 10-55 kA (220 Hz and 230 GHz). Resulting FWHM synthesized beams are listed in Table I. 
Antenna gains were determined from periodic observations of 0528+134 and absolute flux density 
calibration was based on measurements of Uranus. Data were calibrated using the Owens Valley 
software package, MMA, and mapped with AIPS. At all four frequencies, continuum emission is 
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unresolved and peaks at the position of the VLA 43.3 GHz image. Aperture synthesis maps at 89 
and 220 GHz are displayed in Fig. lb and lc. Measured flux densities at all frequencies are listed 
in Table I and displayed in Fig. 2. 



3. Modeling and Discussion 

Our measurements of DO Tau between 8.4 and 230 GHz can be fit by a single power law, 
Fj, oc v a , with index a = 2.39±0.23, shown as a solid line in Fig. 2. Earlier detections of radio 
emission from TTs at wavelengths greater than 1.3 cm yielded values of a between and 1 
(Bieging, Cohen & Schwartz 1984) and have been attributed to free-free emission from ionized 
outflows (Reynolds 1986). Extrapolating from the upper limit of 8.4 GHz emission with a = 1 
(dotted line in figure 2), we estimate that no more than 49% of DO Tau's 43.3 GHz emission 
can arise as free-free radiation from an ionized jet. Fig. 2 suggests the observed mm-wave flux 
originates entirely from circumstellar dust. 

The frequency dependence of the mm-wave dust opacity, (3, can be derived from a, since 
a«2 + /3/(l + A), where A is the ratio of optically thick to optically thin emission from the disk 
(BS, eqn. 1). At frequencies where emission from a circumstellar disk is largely optically thin and 
the Rayleigh-Jeans approximation holds, [3 ~ a — 2. For DO Tau, continuum emission appears to 
be largely optically thin, even at sub-millimeter wavelengths (BS; Mannings & Emerson 1994). 
We estimate 0.39 ± 0.23 for ft 

-7mm) hi good agreement with the BS value of ft). 6— imm> 0-4 ± 0.2. 
There is no evidence for the change in [3 longward of 2mm, postulated by Mundy et al. (1993) for 
a few other TTs. 

An estimate of (3 can also be obtained by fitting the spectral distribution of luminosity 
L u = 4irD 2 vF u with a disk model which takes into account any contribution from optically thick 
emission. Following BSCG and Adams et al. (1990), we assumed power-law radial profiles in disk 
temperature and surface density, T = Tq(R/ Ro)~ q and £ = T,q(R/ Rq)~ p with p = 1.5 or 1.75. 
The millimeter-wave emissivity of the grains, k, is cm g . The outer radius, 

Rd, was allowed to take on values between 22 and 350 AU; the former is the lower limit to disk 
size if all 1 mm emission is optically thick (BS); the latter is the deconvolved half-maximum 
radius of the CO-emitting region from aperture synthesis images. These suggest a disk inclination 
angle, 9, of 40° (Koerner & Sargent 1995). From the 12, 25, and 60 fim IRAS fluxes, which probe 
optically thick regions of the disk, we obtain T = 218 K at 1 AU with q = 0.54, very close to the 
value derived by BSCG for a face-on disk. Best-fit values of (3 and M^, the total disk mass, were 
estimated from the minimum reduced x 2 value. Acceptable fits, with x 2 falling within Ax 2 = 1 of 
its minimum value, were found for our entire range of p and R^ values. The best-fit model, with 
(3 = 0.6 ±0.3, M d = 1.0 ±0.5 x 10~ 2 M , p = 1.75, R d = 350 AU, and x 2 = 0.77, is plotted in Fig. 
3 as a solid line, along with the luminosity distribution derived from IRAS, sub- millimeter, and 
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millimeter observations of DO Tau. Following BSCG (eqn. 20), these parameters yield A « 0.28 
at A = 3 mm and make possible a revised estimate of (3 from the power-law fit to data presented 
here. For A = 0.28 and a = 2.39, we obtain [3 = 0.50 ± 0.23, in good agreement with the value 
obtained from both our disk-model fit and that of Mannings &; Emerson (1994). 

For the ISM, it is commonly assumed that (3 is about 2 in the millimeter wavelength regime 
(Mathis 1990). However, a value of 1.3 has been obtained in recent laboratory studies (Agladze et 
al. 1994). Even lower values are suggested by sub-millimeter observations of T Tauri disks (BS; 
Mannings & Emerson 1994). A variety of explanations have been proposed, including chemical 
composition, physical shape, and grain growth (Wright 1987; BS; Kriigel & Siebenmorgen 1994; 
Ossenkopf & Henning 1994; Pollack et al. 1994). 

For DO Tau, we find (3 « 0.5 and contend that our denser sampling of the sub-millimeter 
regime and spectral coverage extending to longer wavelengths effectively eliminates uncertainties 
that may have complicated other derivations. Grain properties in circumstellar disks are unlikely 
to display the exotic range of chemical composition and physical shapes required to reproduce this 
result in the laboratory. By contrast, the growth of grain size distributions to include particles 
larger than 1 mm accounts for our value of (3 (cf. Miyake and Nakagawa 1993) and is consistent 
with the short theoretical timescales (~ 100 yr) for production of mm-size particles in the early 
solar nebula (cf. Fig. 19, Cuzzi, Dobrovolskis, & Champney 1993). 

If the average grain size in disks steadily increases due to planetesimal formation, [3 should 
decrease monotonically with age. However, DO Tau is relatively young, only a few x 10 5 yrs, 
with an outflow typical of an active disk. Grain growth appears to have already occurred by the 
early T- Tauri phase. Recent 7 mm images of a very young protostar, HH24MMS (Chandler et al. 
1995), also yield a lower value of (3 than found for many older TTs in sub-millimeter surveys (cf. 
BS). These results are inconsistent with a simple picture of gradually decreasing (3; they could 
be explained if mm-size grains grow quickly, followed by generation of a new population of small 
dust grains by planetesimal collisions (Lissauer & Stewart 1993). Long- wavelength observations of 
a statistical sample of TTs disks encompassing a range of ages are clearly required to test this 
hypothesis. 

We are grateful to D. Wood for assistance during the first season of 43 GHz observations at 
the VLA. D.W.K. acknowledges support for this work from NASA grant NGT-51071. The Owens 
Valley millimeter-wave array is supported by NSF grant AST-9314079. Research by A.I.S. on 
protoplanetary disks is furthered by NASA grant NAGW-4030 from the "Origins of Solar Systems" 
program. The research described in this paper was carried out by the Jet Propulsion Laboratory, 
California Institute of Technology, and was sponsored by a fellowship from the National Research 
Council and the National Aeronautics and Space Administration. 
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Frequency 


Flux Density 


Statistical 


Total 


Synthesized 


Beam 


Parameters 


(GHz) 


(mJy) 


error (mJy) 


error (mJy) 


Pmaj 


Bmin 


PA 


8.4 


<0.17 (3cr) 






0.41" 


0.38" 


120° 


22.5 


<0.76 (3cr) 






0.25" 


0.23" 


-10° 


43.3 


1.80 


0.35 


0.71 


0.68" 


0.53" 


-78° 


89.2 


14.2 


0.9 


3.74 


2.86" 


2.10" 


82° 


111.2 


30.2 


4.5 


10.5 


13.0" 


5.4" 


-69° 


221.5 


98.6 


4.9 


24.2 


3.98" 


3.15" 


62° 


232.0 


137.5 


4.9 


32.4 


3.41" 


3.15" 


-84° 



Table 1: Radio (VLA) and mm-wave (OVRO) Continuum Flux Densities from DO Tauri. Total 
errors include 20% uncertainty in the absolute flux density calibration. 
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